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ABSTRACT 

We present the aim and first results of the RApid Temporal Survey (RATS) made 
using the Wide Field Camera on the Isaac Newton Telescope. Our initial survey covers 
3 square degrees, reaches a depth of V~22.5 and is sensitive to variations on timescales 
as short as 2 minutes: this is a new parameter space. Each field was observed for over 
2 hours in white light, with 12 fields being observed in total. Our initial analysis 
finds 46 targets which show significant variations. Around half of these systems show 
quasi-sinusoidal variations: we believe they are contact or short period binaries. We 
find 4 systems which show variations on a timescale less than 1 hour. The shortest 
period system has a period of 374 sec. We find two systems which show a total eclipse. 
Further photometric and spectroscopic observations are required to fully identify the 
nature of these systems. We outline our future plans and objectives. 

Key words: methods: data analysis - techniques: photometric - surveys - stars: 
variables: - stars: binaries - general - Galaxy: stellar content 


1 INTRODUCTION 

The intensity of stellar objects can vary on a wide range 
of time-scales, ranging from seconds to months to years. A 
large number of projects now exist whose aim is to detect 
such varying sources. The reasons for this are many, but 
include the search for extra-solar planets and interacting 
binary stars. Most of these surveys are sensitive to timescales 
longer than a day. It is only recently that such surveys have 
been sensitive to shorter term timescales. For instance, the 
0.3° x0.3° survey using the Canada France Hawaii Telescope 
was sensitive to variations on timescales as short as ~8 min 
(Lott et al 2002). While the Faint Sky Variability Survey 
has a much larger survey area, it was sensitive to variations 
only as short as ~24 min (Groot et al 2003, Morales-Rueda 
et al 2004). In principle the Super Wasp project is sensitive 
to variations on timescales as short as a few mins (Christian 
et al 2004). However, they are sensitive to relatively bright 
objects, V ~7-15. 

Why is it important that we extend the parameter 
search down to periods shorter than 10 min? Recently a 
new class of object has been discovered in which coherent 
intensity variations have been detected on timescales of ~10 
min or less, with the shortest being 5.4 mins (see Cropper 
et al 2004 for a review). It is thought that these systems are 
interacting white dwarf-white dwarf pairs which have no ac¬ 
cretion disc, and the observed period represents the binary 


orbital period. As such, they are expected to be amongst the 
first sources to be detected using LISA, the planned gravita¬ 
tional wave observatory (Nelemans, Yungleson & Portegies 
Zwart 2004). 

These systems are at the short period end of the pe¬ 
riod distribution of white dwarf-white dwarf binaries, or 
AM CVn systems. For orbital periods less than 80 min, the 
secondary (mass donating) star cannot be a main sequence 
star. Further, for periods shorter than 30 min, the secondary 
must be a Helium white dwarf (Rappaport, Joss & Webbink 
1982). There are around 13 of these AM CVn systems cur¬ 
rently known. 

Theoretical predictions suggest that AM CVn systems 
are commonplace (eg Nelemans et al 2001). It is not clear 
if the discrepancy between the observed number and pre¬ 
dicted number of systems is due to an over estimate in the 
predicted number or if many more systems await discovery. 
To help resolve this question we have started a project, the 
RApid Temporal Survey (RATS), whose aim is to discover 
new objects which vary in a coherent manner on timescales 
less than ~lhr. As a significant by-product of this search we 
expect to discover many new objects which vary on longer 
periods. This paper presents our strategy for discovering new 
AM CVn systems; the positions and light curves of those 
systems which we have discovered as part of our initial set 
of observations; a discussion regarding the nature of these 
sources and a brief discussion as to what these observations 
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Table 1. The field centers for our observations. The co-ordinates 
are in J2000 

imply regarding the population density of AM CVn systems 
and other objects. 


2 OBSERVATIONS AND REDUCTION 

2.1 Observational strategy 

Our observations were made over 3 nights starting on 28 Nov 

2003 using the 2.5m Isaac Newton Telescope on La Palma 

and the Wide Field Camera (WFC). The WFC consists of 

// 

4 EEV 2kx4k CCDs, with each pixel corresponding to 0.33 
on the sky. The approximate field of view is 33 x 33 (with 
a 11 x 11 gap missing on one of the corners). The orien¬ 
tation of the instrument was such that our field center was 
located near the center of CCD 4 (see Groot et al 2003 for 
a schematic diagram showing the orientation of the chips). 

Since our goal was to obtain photometry with the high¬ 
est possible time resolution, we did not use a filter - eg they 
were ‘white light’ observations: since our observations were 
made near the zenith, differential atmospheric diffraction 
was reduced to a minimum. However, at the start of each 
observation we also obtained BVI images. Since we did not 
take observations of standard stars before every field obser¬ 
vation, the BVI magnitudes reported in Table 2 are prob¬ 
ably accurate to only ~0.1 mag. Our exposure times were 
30 sec in duration, with another ~42 sec of readout time for 
the whole array. Each field was observed for around 2-2.5 
hrs. We obtained bias frames and flat fields on each night. 

2.2 Field selection 

Our aim was to select fields which were close to the galactic 
plane so that there were as many stars in the field as pos¬ 
sible, but not so close to the plane that crowding became 
a serious issue. Further, we wanted to select fields that did 
not contain stars brighter than ~10 mag, which might cause 
charge overflow on the CCDs. To make this selection we ob¬ 
tained the APM catalogue 1 for a range in galactic latitudes 
(5 < |6| < 25) and wrote a routine which optimizes the field 

1 http://www.ast.cam.ac.uk/~apmcat 


selections for the given sky coordinate range and the shape 
and size of the field of view. For an observing run lasting n 
nights, our optimization routine finds the same number (n) 
of fields around 4 positions in the sky that pass close to the 
zenith during the night with about 2 hr spacing. 

This is done using an algorithm that first places n ran¬ 
dom fields in the sky in a 5x5 degree box centered around 
a given pointing. Then a simulated annealing based opti¬ 
mization routine is used to find a combination of n non¬ 
overlapping pointing’s within this 5x5 degree box. The op¬ 
timization tries simultaneously to maximize the number of 
V= 15-20 mag stars and minimize the number of brighter 
stars included. There is also a hard limit, so that no stars 
brighter than V=10 mag should appear in any of the fields. 
So, as a result we will get, for each night, 4 fields per night 
that are observable for around 2hr close to the zenith after 
each other, thus filling all the nights optimally. 

2.3 Data Reduction 

A series of ~100-120 30 sec exposures were made of each 
field: auto-guiding was used so there is no large x, y shift 
between frames. For each chip we flat-fielded and bias sub¬ 
tracted each image. The analysis (for each chip) consisted 
of following steps: the first 5 frames of each field were com¬ 
bined in order to form a ‘master’ image of the field. This 
master image was then searched for all the sources above 
a 5 a detection limit. Next, we determined for each of the 
remaining frames an accurate image shift in relation to the 
master image. This was done in order to be able to use small 
enough apertures for light curve extraction, thus maximiz¬ 
ing the S/N ratio for the photometry. In practice, this was 
achieved by using the positions of a couple of bright (but not 
saturated) stars. After having determined (and applied) the 
shifts for individual images, we obtained quasi-differential 
photometry using the median brightness of the sources in 
each frame as a comparison (to remove any variations in the 
sky transparency) - this was done for all the images of the 
field. In total around 33000 stars were detected in our survey 
which covered a total of 3 square degrees. 

All the resulting light curves were analysed and those 
which had bad magnitude values or errors above 0.1 mag or 
non-detections in any frame, were discarded. This reduced 
the total number of light curves by around 2000. All the 
‘good’ light curves were analysed using a Lomb-Scargle pe¬ 
riod search algorithm, and those light curves which showed 
a significant peak in their power spectrum were flagged for 
closer inspection. The significance of the highest peak in 
each power spectrum was estimated in the following man¬ 
ner. First, we took the highest peak power value and mean 
power value for each power spectrum. Then we defined the 
significance statistic of the highest peak in the power spec¬ 
trum by computing the ratio of it to the mean power value. 
The final selection of ‘interesting’ variable targets was car¬ 
ried out manually. For each chip we produced a plot of ap¬ 
parent magnitude against the significance statistic for all the 
stars. When the overall number of fields is not too large (as 
in the case here), such a plot provides a convenient way of 
of identifying the likely variables in each field. 

To test our procedure we selected the field of RX 
J0806+15, the 5.4 min ultra-compact binary system. Its 
brightness variation is highly stable with a peak-to-peak 
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Figure 1. The results of our routine to detect sources varying 
in their intensity using the field of the ultra-compact binary RX 
J0806+15 as a test case. The magnitudes have been scaled so 
they are approximately that of V mag and higher significance 
values imply strong variability. Our routine clearly identifies RX 
J0806+15 as being highly variable. 

modulation of ~0.3 mag in white light (Ramsay, Hakala & 
Cropper 2002). We ran our data analysis procedure in the 
manner outlined above. For each star in the field we plot 
the significance of the highest peak in its power spectra (as 
defined above) as a function of brightness in Figure Q RX 
J0806+15 was easily picked up by our routine. We are there¬ 
fore confident that we can detect sources which have coher¬ 
ent variability with amplitudes of ^10 percent or lower. 

For those sources identified as being significantly vari¬ 
able, we used autophotom (Eaton, Draper & Allan 2003) to 
obtain ‘optimal’ photometry for each source together with 
suitable comparison sources. Whilst our procedure for iden¬ 
tifying variable sources is relatively simple, we have demon¬ 
strated that it will detect targets such as RX J0806+15 with 
ease. However, it is probable that additional variable sources 
still await detection using a more sophisticated detection 
routine. 

2.4 Source positions 

To determine accurate positions for our variable sources, we 
compared our fields with the Digital Sky Survey images and 
the Hubble Guide Star Catalogue. We typically used 8 refer¬ 
ence stars for each chip. We then compared the x, y position 
of each reference stars in our images and used astrom (Wal¬ 
lace & Gray 2002) to obtain J2000 co-ordinates for each of 
our sources.The residuals to the fits were typically ~ 0.5 . 
We searched for known sources in the SIMBAD database 
using a aperture of radius 10 centered on our source po¬ 
sitions: there were no known sources at these positions. All 
our sources are therefore newly discovered variable sources. 


3 INITIAL RESULTS 

Out of the 33000 stellar objects in our survey, we have identi¬ 
fied 46 which have showed evidence for significant variability. 
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Figure 2. Top panel: The number of variable source discovered 
in our survey as a function of magnitude. Bottom panel: The 
percentage of variable sources as a function of the number of 
stars in a magnitude range. 

Their position, V mag and colour and the characteristics of 
their light curves are shown in Table [3 

We show in Figure[21the actual number of sources which 
varied as a function of apparent brightness and also the frac¬ 
tion of sources which were found to be variable as a func¬ 
tion of brightness. We find variable sources in the range 
15 < V < 22, with an approximately uniform spread in 
the number of sources as a function of apparent brightness. 
Although, the number of sources in the bins are relatively 
small, we note that the fraction of sources which are variable 
is not uniform, showing a peak in the bin 15.5 < V < 16.5. 
We discuss one possible reason for this in §4. 

The light curves of our 46 variable sources are shown in 
Figures 3-5. They can be split up into 5 broad categories: 
those showing sinusoidal-like variations; coherent variations; 
eclipse-like behaviour; flaring; and more irregular behaviour. 
We now go onto discuss various objects in more detail. 

3.1 Objects varying on periods of less than lhour 

The main goal of this project is to discover sources which 
vary in their intensity in a coherent manner on timescales of 
~1 hour or less. We have discovered one source with a coher¬ 
ent period of 374 sec, 2 systems showing sinusoidal variations 
on timescales of ~1 hour, one which shows a variation on 40 
or 80 mins and one with a possible variation on a period of 
75 mins. 

The system with the shortest period was RAT 
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Figure 3. Light curves of all interesting sources. 
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Figure 5. Light curves of all interesting sources (cont). 


Folded on 374.2 sec 



Figure 6. White light data of RAT J0455+1305 folded on the 
best fit period of 374 sec. 


J0455+1305 which was found to vary on a period of 374 sec 
and an amplitude of 0.13 mag. The source is detected in the 
USNO-B1 Catalogue, but not in the ROSAT all-sky survey. 
The folded light curve appears sinusoidal (Figure |SJ|. The 
shape and amplitude are virtually identical to the discless 
ultra-compact binaries mentioned in §1. On the other hand 
pulsating white dwarfs and sub-dwarf pulsators (EC 14026 
stars) also show such coherent modulations. Dedicated ob¬ 
servations of this object will be presented in a future paper. 

The two systems which show sinusoidal variations 
on periods close to lhr, RAT J0455+1254 and RAT 
J0807+1507, have amplitudes of 0.06 and 0.25 mag respec¬ 
tively. The latter has a light curve which shows clearly 
peaked maxima. If their periods are indeed close to 1 hr, 
then this would rule out them being W UMa systems - con¬ 
tact binaries which have orbital periods in the range 0.2-1.5 
days. 

RAT J0449+1756 shows clear variations with an ampli¬ 
tude of ~0.04 mag. It is unclear whether there is a modula¬ 
tion on a period of 40 or 80 min. Unfortunately this object 
was just off the chip edge when we made the BVI images. 
It was, however, detected in the 2MASS survey: it shows 


JHK colours which are blue. Its blue colour and its period 
are similar to that of the AM CVn systems: in this case, it 
would be expected to show strong Helium lines in emission 
in its optical spectrum. 


3.2 Objects showing eclipse-like features 

There are three sources which show eclipse-like behaviour: 
RAT J0740+2340 (a total eclipse of depth 0.1 mag), RAT 
J2310+3414 (an eclipse-like feature of depth ++0.2 mag) and 
RAT J2257+3424 (a partial eclipse with depth 0.2 mag). 
Each system shows red optical colours out of eclipse (cf Ta¬ 
ble 2) and were all detected in the 2MASS survey 2 : RAT 
J0740+2340 has colours of a M4-5 dwarf star while RAT 
J2257+3424 is consistent with an early K dwarf star. In the 
case of RAT J0740+2340 the white light and BVI images 
show a slight elongation of the source implying the colours 
for this source may be contaminated by a close ‘companion’ 
star. RAT J2310+3414 shows infrared colours, (H-K)=0.52, 
(J-H)=-0.08, which are not consistent with that of main se¬ 
quence stars. 

The eclipse of RAT J0740+2340 takes ~20 min to de¬ 
scend into and out of eclipse. The total duration of the 
eclipse is ~65 min (half eclipse depth). The rather low eclipse 
depth rules out it being a cataclysmic variable which would 
show a much greater eclipse depth. RAT J2257+3424 maybe 
an Algol system in which we have detected the secondary 
eclipse. 

As the shape of the eclipse profile of RAT J0740+2340 
closely resembles one produced by a planetary transit, we 
have studied it’s properties further. Assuming the eclipse 
is produced by a transiting planet, we can work out some 
preliminary parameters for the system. Firstly, the 0.1m 
depth of the eclipse implies that the radius of the planet 


2 http://pegasus.phast.umass.edu 
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Table 2. Those sources discovered in the first reduction of our INT/WFC survey data. We quote the period of any source showing 
periodic behaviour, together with the change in white light magnitude and the main characteristic of the light curve (sin: sinusodial 
in shape; eel: eclipse; shortp: short period variability). *RAT J0208+3605 is ~ 1 distant from another source from which it was not 
possible to fully resolve: the quoted magnitudes are the combined. 


has to be ~0.3 of the stellar radius. Secondly, the ratio of 
ingress/egress lengths to the duration of the total part of 
the eclipse implies that the transit occurs at an apparent 
stellar latitude of about 30°. 

Some preliminary spectra of RAT J0740+2340 show 
that the star is of very late spectral type. Assuming a M4- 
6V spectral type, this would imply a stellar radius of ~0.3 
Rq and planetary radius of about 70000 km, i.e. the same 
as Jupiter. Using the stellar radius and the implied transit 
latitude together with the FWHM eclipse duration yields 


an orbital velocity of 90 km/sec, which in turn leads to an 
orbit of roughly 3 days. This is similar to that observed in a 
handful of systems (eg Mazeh, Zucker & Pont 2005). How¬ 
ever, further observations are required to confirm the orbital 
period of this object. 


3.3 Objects showing longer term variations 

Approximately half of our newly discovered sources show si¬ 
nusoidal or quasi-sinusoidal intensity variations with periods 
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greater than ~lhr. Sources which exhibit these characteris¬ 
tics include pulsating stars such as <5 Scuti stars, and the W 
UMa and Algol binary systems. 

<5 Scuti stars show intensity variations on timescales be¬ 
tween 30 min and ~8 hrs, with amplitude variations less 
than ~1 mag. However, these systems tend to show asym¬ 
metric light curves, with the rise to maximum taking a 
shorter time than the descent from maximum. We have 
discovered one candidate 5 Scuti system: RAT J0305-0036 
shows an asymmetric light curve; a period of ^2 hrs and an 
amplitude of ~0.3 mag. 

W UMa systems are well known systems, being thought 
to account for 1 in every 500 stars (Rucinski 2002) and 
show orbital periods between 4hrs - 1.5 days. Other sys¬ 
tems include binary systems in which the mass losing star 
is irradiated by the mass accreting star: this results in the 
companion star being apparently brighter when it is viewed 
face-on compared to when we observe it from the rear. In 
recent years it has become clear that many sub-dwarf B 
stars (sdB) are in such binary systems: these stars are Ex¬ 
treme Horizontal Branch stars which have Helium cores and 
only a thin layer of Hydrogen on their surface (eg Maxted 
et al 2002). The companion star in these systems could be 
main sequence stars, white dwarfs or even brown dwarf stars. 
Those systems with white dwarfs have potentially wide as- 
trophysically significance since they maybe SN la progeni¬ 
tors (Maxted, Marsh & North 2000). We expect that many 
of our newly discovered sources will be contact and non- 
contact binaries. 


3.4 Flare sources 

We detected two flare like objects in our survey. RAT 
J0456+1254 showed a 0.25 mag flare while RAT J0806+1540 
showed a ~1 mag flare. RAT J0456 was detected in all three 
BVI bands while RAT J0806 was detected only in the / 
band. The former has colours consistent with a main se¬ 
quence M5/M6 star (Pickles & van der Kruit 1990), which 
using the absolute magnitude relation in Allen (1973) gives 
a distance of ~725 pc for a M5V star and ~380pc for a M6V 
star. RAT J0806+1540 is therefore likely to be an even later 
type dwarf star or more distant than RAT J0456. M dwarf 
stars are known to show prominent flare behaviour. 

3.5 Minor Planets 

As expected a number of minor planets were detected: a 
total of 16 were identified. Their positions were determined 
for 3 epochs using the procedure in D these were passed 
to our minor planet colleagues based at the University of 
Helsinki. They determined if there were any known minor 
planets at positions at these epochs - 5 were not previously 
known. 

To accurately determine the orbit of a minor planet, and 
therefore give rise to a designation for that system, positions 
have to be determined on several nights. Unfortunately, since 
the systems were not detected until after the observing run 
was over, we did not get these positions. For future surveys 
we will run our source detection software for each field in 
near real-time. This will allow any new minor planet to be 
followed up on a subsequent night. 


4 DISCUSSION & SUMMARY 

We have presented the first results of our RApid Temporal 
Survey which covered 12 INT/WFC fields, giving a total sky 
coverage of approximately 3 square degrees. Our survey ex¬ 
plores a new parameter space - namely searching for objects 
which show intensity variations on timescales as short as a 
few mins. Using our relatively simple algorithm to detect 
such sources we have discovered 46 sources which showed 
variability on timescales between 374 sec and 4-5 hrs. For 
sources which vary on periods longer than ~8 hrs, it is prob¬ 
able that we would not have identified them as such in our 
short observing window. To determine the nature of the new 
variable sources, plans are in place to obtain followup pho¬ 
tometry and spectroscopy. 

We estimate that our white light images reached a depth 
corresponding to an equivalent V~22.5 (the exact depth de¬ 
pends on the colour of the source). The distance to the edge 
of the galaxy on the far side of the galactic center in the 
plane is ~22.5kpc. For a galactic latitude of 20° (typical 
of fields observed in our initial survey) we can detect every 
star with a spectral type earlier than ~K2 out to the edge 
of the galaxy. However, the scale height of the galaxy is of 
the order of 200pc. This implies that for stars within 3 scale 
heights of the galaxy we can detect every star as faint as 
~M5. 

We noted in J]]tliat there was an peak in the percentage 
of objects which were variable in the range 15.5< V < 16.5. 
If we assume that at a distance of 3 scale heights the num¬ 
ber of additional sources which we detect is negligible then 
that implies a distance of 1.7kpc for galactic latitude of 20°. 
For U=16 and a distance of 1.7kpc this implies a depth of 
My= 6. We note that W UMa systems have a range in ab¬ 
solute magnitude of My =2-6 with systems with the short¬ 
est orbital periods (0.2 days) being intrinsically fainter than 
longer period systems. We suggest that if the peak in the 
distribution of percentage of objects which are variable is 
real (to be confirmed in the extension of the survey) then it 
is due to W UMa systems: after V ~ 16 there are no more 
of these systems to discover. This is consistent with the sug¬ 
gestion that a significant fraction of our newly discovered 
systems are W UMa systems. 

One of the main science goals of this project is to search 
for AM CVn systems. Although this paper presents only our 
initial findings, we have found only one or two (at most) 
systems which could, with further observations, be classed 
as candidate AM CVn systems. How many were expected 
from population synthesis models? Nelemans et al (2001) 
predicted a space density of 2.1xl0 -4 pc -3 which was re¬ 
vised down-wards by around 60 percent by Nelemans, Yun- 
gelson & Portegies Zwart (2004). The distance that we can 
detect such systems is dependent on their absolute mag¬ 
nitude which in turn is dependent on their orbital period. 
For longer period systems such as GP Com, My ^12, im¬ 
plying we could detect them out to a distance of 1.6 kpc 
(we expect to detect systems with shorter orbital periods 
out to greater distance). If we take a uniform distribution 
of AM CVn systems, then in a 3 square degree survey the 
above space density predictions imply around 150 systems 
(with an uncertainty of at least 2). However, Nelemans et 
al (2004) made an investigation regarding selection effects 
and found that for an optical limit of V = 20 only approxi- 
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mately 1 percent of systems with orbital periods less than 25 
mins would be detected as an AM CVn system in the opti¬ 
cal band, implying only 1 -2 systems in our 3 square square 
degree survey. Since our survey extends to a fainter opti¬ 
cal limit and work is needed to compare the overall space 
density at different orbital periods, it is too early to make 
definitive conclusions regarding the comparison of theoreti¬ 
cal predictions and observed space density estimates. 

In conclusion, we have found by exploring a new pa¬ 
rameter space around 46 new objects which have not before 
been identified as variable sources. How does this discovery 
rate compare with other surveys? Mukadam et al (2004), 
for instance, found 35 pulsating white dwarfs in 125 nights 
of observing time. We have discovered one pulsating star in 
3 nights, giving a similar discovery rate. However, we have 
found many new objects in addition to this one very short 
period system. More observational data is needed to con¬ 
firm their nature: this work is in hand, and our results will 
be presented in a followup paper. 
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